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The much debated issue of the transverse single spin asymmetry An observed in the inclusive 
large Pt production of a single hadron inpp interactions, p^ p — > irX, is considered in a TMD 
factorization scheme. A previous result 1, 2\ stating that the maximum contribution of the Collins 
effect is strongly suppressed, is revisited, correcting a numerical error. New estimates are given, 
adopting the Collins functions recently extracted from SIDIS and e + e~ data, and phenomenological 
consequences are discussed. 

PACS numbers: 13.88. +e, 12.38.Bx, 13.85.Ni 



I. INTRODUCTION AND FORMALISM 



The understanding of spin effects, and in particular, for parity conserving processes, of transverse Single Spin 
Asymmetries (SSAs) has always been one of the major challenges for QCD and any fundamental quantum field 
theory. Such effects, abundantly observed in experiments, are not generated by the perturbative dynamical properties 
of the Standard Model elementary interactions and originate from more profound phenomena related to the intrinsic 
structure of the nucleons. Thus, the study of these often neglected effects has recently opened a new phase in our 
exploration of the partonic composition of hadrons. 

Most progress has occurred in the study of the rich azimuthal dependences measured in Semi-Inclusive Deep Inelastic 
Scattering (SIDIS) of leptons off polarized nucleons by dedicated experiments, HERMES (DESY), COMPASS (CERN) 
and JLab. These SIDIS azimuthal asymmetries arc interpreted and discussed in terms of new, unintcgratcd, Transverse 
Momentum Dependent distribution and fragmentation functions (shortly, TMDs); these offer new information on the 
properties of quarks and gluons, which go well beyond the usual one-dimensional description of Partonic Distribution 
Functions (PDFs) in terms of longitudinal momentum fraction only. 

In particular the Sivers distributions 043 and the Collins fragmentation functions @ have been extracted 043 
from SIDIS data, and, thanks to complementary information from Belle on the Collins function [l]| [l4|], a first 
extraction of the transversity distribution has been possible (l5l fl(| . 

All these analyses have been performed in the 7* — p center of mass (cm.) frame, within a QCD factorization 
scheme, according to which the SIDIS cross section is written as a convolution of TMDs and elementary interactions: 

da e P ^'hx = J2f\ /p (x,k ±] Q 2 )®da i ^®D h/q (z, P± -,Q 2 ), (1) 

where k± and are, respectively, the transverse momentum of the quark in the proton and of the final hadron with 
respect to the fragmenting quark. At order k±/Q the observed Pt of the hadron is given by 

P T = zk ± +p ± . (2) 

There is a general consensus fl7l - fl9j that such a scheme holds in the kinematical region defined by 

P T ~ k± =i Aqcd < Q ■ (3) 

The presence of the two scales, small Pt and large Q, allows to identify the contribution from the unintegrated partonic 
distribution (Pt — k±), while remaining in the region of validity of the QCD parton model. At larger values of Pt 
other mechanisms, like quark-gluon correlations and higher order perturbative QCD contributions become important 
[L9j - |2ll ]. A similar situation [l8l. [22| - [27| holds for Drell-Yan processes, AB — > £ + £~X, where the two scales arc the 
small transverse momentum, qt, and the large invariant mass, M, of the dilepton pair. 
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The situation is not so clear for processes in which only one large scale is detected, like the inclusive production, at 
large Br, of a single particle in hadronic interactions, AB — > CX. However, the most striking and large SSAs have 
been [28W31I ] and keep being measured [3214371 ] in these reactions. The TMD factorization for these processes was first 
suggested in Refs. [3 l ^ and adopted in Refs. [38l - |4(j| to explain the large single spin asymmetries observed by the 
E704 Collaboration 2{|. The same approach led to successful predictions [4l|, l42[ for the values of An measured at 
RHIC 0. 

Alternative approaches to explain the origin of SSAs, linking collincar partonic dynamics to higher-twist quark- 
gluon correlations, were originally proposed in Refs. [44| - l48| and phenomcnologically used in Refs. |49l- [ 52 [l . These 
two approaches, the TMD factorization and the higher-twist correlations, have been shown to be related [53|, HH and 
consistent with each other [23|, [24|, Hf| . 

However, a definite proof of the validity of the TMD factorization for hadronic inclusive processes with one large 
scale only is still lacking. Due to this, the study of dijet production at large Bt in hadronic processes was proposed 
[561 - 459( 1 . where the second small scale is the total qr of the two jets, which is of the order of the intrinsic partonic 
momentum This approach leads to a modified TMD factorization approach, with the inclusion in the elementary 
processes of gauge link color factors [60l - l62l |. Despite the identification of two separate scales, some problems with 
the TMD factorization for hadronic processes with the final inclusive production of two jets, or two hadrons, have 
been recently pointed out [63|, |64j]. TMD factorization is expected to work for the observation, inside a jet with 
large transverse momentum Br, of a final hadron with a transverse momentum with respect to the jet direction, like 
p^p -> jet + 7T + X [65l-l67t. 

In this paper we consider SSAs in p'p — > ir X processes, with only one large Bt final pion detected, for which data 
are available. We adopt the TMD factorization scheme 



ss, witn only 



do**^"* = ^ fa/p(xa,k^ a ;Q')^f b/p (x b ,k ±b ;Q')^da ab ^ cd ^D n/c (z,p ±] Q'), (4) 

a,b,c,d 

as a natural phcnomenological extension of the corresponding collincar factorization, based on the convolution of 
integrated parton distributions (PDFs) and fragmentation functions (FFs) with QCD elementary dynamics; this 
collincar factorization works well in computing unpolarized cross sections, but fails in explaining SSAs as there is no 
single spin effect in the collinear PDFs and FFs, and in lowest order QCD dynamics. In the polarized case, for a 
generic process (A, Sa) + (B, Sb) — > C + X, Eq. (@| explicitely reads Q: 



E c da( A ' S ^+( B ' S ^ c+x _ yv f dx a dx b dz 

a,b : c,d,{\} 



d 3 P c 



—\ 5- d 2 k ±a d 2 k ±b d 3 p ± 5(p ± - p c ) J{p±) 

IblT^XaXbZ^S 



P\{ A \f A fa/A, 3a ( x a, k± a ) Px^* fb/B,S B ( X b, k ± b ) (5) 

is * 



Further details and a full explanation of the notations can be found in Refs. [H, (where p ± is denoted as fej_cj 

There are many contributions, from different TMDs, to the SSAs according to the above expression. In Rcf. [1| it 
was argued that only the Sivers effect contributes significantly; a further small contribution from the Collins effect is 
possible (largely suppressed by phase integrations), while all other TMD contributions are utterly negligible. However 
that conclusion was affected by a wrong sign jga, |69| in one of the elementary interactions and the Collins effect 
contribution was underestimated. It remains true that all other contributions to the SSAs are negligible. 

Reconsideration of the Collins contribution 

We reconsider here the Collins contribution to the SSA 

A n = tttti where dcr 1 ^ = — • (6) 

da 1 + da* d 6 p^ 

Such a contribution can be computed according to the TMD factorized expression [l|, Q : 

[da' - da+] CoUiDS = V / dx * dx » dz d 2 k±a d 2 k±b d 3 p± §{p± . p ) J(p±) S{s + 

^ ,J 16 TT Z X a X b Z Z S 

q a , b ,q c ,d 

x A T q a (x a , k± a ) cos((j> a +ipi- ip 2 + <f>%) (7) 
x f b/p (x b ,k ±b ) \M<tM$] A N D ,(z, P± ), 
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which can be easily interpreted: the transversity distribution Axq a (or h\ a ) of the quark q a , couples with the unpolar- 
izcd TMD for parton b, f b / p , and the Collins fragmentation function A N D iK , q t = (2p±/zMh) H^ qc of quark q c into 
a 7T. The M® are the independent helicity amplitudes defined in Refs. [HQ, describing the lowest order QCD inter- 
actions and the quantity M® is proportional to the elementary spin transfer cross-section da q 



Ib-t-qld 



da q 



The (suppressing) phase factor cos(0 a + <p\ — ip2 + 4>%) originates from the k± dependence of the unintegrated 
transversity distribution, the polarized elementary interaction and the spin-pj_ correlation in the Collins function. 
The explicit expressions of ipi , ip 2 and <p^ in terms of the integration variables can be found via Eqs. (60)- (63) in 
and Eqs. (35)-(42) in Q. 

The denominator of Eq. ([6]) is twice the unpolarized cross-section and is given in our TMD factorization scheme by: 



[<!** + do*] = 2da unp = 2 



da pp ~ 



■x 



- E 



d 3 P« 

dx a dx b dz 

16 7T 2 X a Xb Z 2 S 



d 2 kj_ a d 2 kj_ b d 3 p ± S(p ± -p c ) J(p x ) S(s + i+u) 



(8) 



x f a/p (x a> k±a) f b/p (x b , k ±b ) |M°| 2 + |M 2 °| 2 + |M 3 °| 2 



ab—tcd 



The explicit expressions of Mj M®, which give the QCD dynamics in Eq. (0, can be found, for all possible elementary 
interactions, in Ref. Q. Unfortunately, it turns out that, for the important qg — > qg channel, as given in the last line 
of Eq. (55) in Ref. [l| and in the last term of Eq. (71) in Ref. Q, there is an overall wrong sign. Moreover, at the time 
of the numerical estimates of Ref. [l[ , no information was available on the transversity distribution and the Collins 
functions, and we computed the maximum Collins contribution to An by replacing these unknown functions with their 
upper limits imposed by positivity conditions. Much more information is now available on the transversity distribution 
and the Collins fragmentation function from studies of SIDIS data by COMPASS and HERMES Collaborations and 
e + e~ data by Belle Collaboration [H|,[l6j]. Notice that the Collins fragmentation function is expected to be process 
independent [541 . l66l [7l| , so that this new information can be used in the inclusive hadronic reactions we consider 
here. 

Due to these reasons - the correction of a previous numerical error, and the use of actual information on the 
functions contributing to Eq. - we are now able to attempt new realistic estimates of the contribution of the 
Collins effect to the intriguing data on Am- 



II. THE ROLE OF THE COLLINS EFFECT 



Thus, in this section, we investigate in much detail the phenomenology of the Collins effect for pion SSAs at STAR, 
BRAHMS and E704 kinematics, according to Eqs. ©-(H]), exploiting all the available and updated information on 
the transversity distribution and the Collins fragmentation functions coming from SIDIS and e + e~ annihilation data. 
In addition, we correct the numerical error in one of the elementary interactions mentioned above, which affected the 
conclusions of Ref. pj . 

The first combined extraction of the quark transversity distribution and the Collins function was presented in 



Ref. [15j. We will refer to it as the "SIDIS-1" fit. In this extraction, the Kretzer set of unpolarized FFs [72( was 
adopted. An updated extraction of the transversity and Collins functions was presented in Ref. [l6|. We will refer 
to the corresponding set of parameterizations as the "SIDIS-2" fit. In this case, the set of pion and kaon FFs by de 
Florian, Sassot and Stratmann [73j, which became available at that time, was considered. 

Let us recall the main features of the parameterizations adopted in Refs. [IHOjl. The analysis of SIDIS and e + e~ 
data is performed at leading order, 0(k±/Q), in the TMD factorization approach, where Q is the large scale in the 
process. A simple factorized form of the TMD functions was adopted, using a Gaussian shape for their kj_ dependent 
component. For the unpolarized parton distribution and fragmentation functions we have: 

e -k 2 ± /{k 2 ± ) e -P 2 ±/(pl) 

f q/p (x,k±) =f q / p (x) 2 , D h/q (z,p±) = D h/q {z) 2 , (9) 

n \ K x) n xP±) 

where (fej_) and (pj_) have been fixed by analyzing the Cahn effect in unpolarized SIDIS processes, see Ref. 0: 

(k 2 x ) = 0.25 GeV 2 , (pi) = 0.20 GeV 2 . (10) 



For the usual integrated PDFs f„/ v (x) we adopted the GRV98 set [74j and, as said above, for the integrated FFs 
D h/q (z) we used the Kretzer set [72| for the SIDIS-1 fit and the DSS one [73| for the SIDIS-2 fit. We have taken into 
account their DGLAP QCD evolution. 
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The quark transversity distribution, Axq(x, k±), and the Collins fragmentation function, A N D h / q t(z,p±), have 
been parametrized as follows: 

A T q(x,k ± ) = -N q T (x) [f q/p (x) + Aq(x)] G , (11) 

A N D h/q ,(z, P± ) = 2Af q c (z)D h/q (z)h(p ± ) T~2\ — - (12) 

where Aq(x) is the usual collinear quark helicity distribution, 



and 



with \N q {C} \ < 1. Moreover, 



Af q c (z) = iVf - *)'. ^L2^lL s , (14) 

7? V 



h(p 1 _) = V2-e^-e-A/ M Z. (15) 

With these choices, the transversity and Collins functions automatically fulfill their proper Soffer and positivity 
bounds respectively, for any values of the (x,k±) and (z,p±) variables. The quark helicity distributions Aq(x), 
required for the Soffer bound, are taken from Ref. [75|. The term [f q / p (x) + Aq(x)] in Eq. (jlip is evaluated at the 
initial scale and evolved at the appropriate Q 2 values using the transversity evolution kernel. Similarly, for the Q 2 
evolution of the Collins function, which remains so far unknown, we considered the unpolarized DGLAP evolution of 
its collinear factor D h / q (z). 

Despite the simplicity of these functional forms, they still involve, in the most general case, a huge number of free 
parameters. In Refs. (l5l . [iH ] we therefore adopted some additional, physically motivated assumptions, in order to 
keep the number of free parameters reasonably low. First of all, for the transversity distribution we used only valence 
quark contributions. In addition, for the fragmentation functions we considered two different expressions for M q , 
corresponding to the so-called "favoured" and "unfavoured" FFs, J\f^ v (z) and Af^ nf (z); for example, for pions, we had: 

Af? +/u (z)=M° +/3 (z)=AfZ- /u (z)=M?- /d (z)=Mg(z), (16) 



^ +/a (z) = N^ /d (z) =AC C /u (z) =N?_ /d (z) =N? ±/a (z) = NZ±,M (17) 

Notice, however, that our complete parameterization of the Collins FFs, Eq. (fT2"|. allows for further differences 
among parton flavours, possibly contained in the usual unpolarized FFs. 

In addition, we kept a flavour dependence in the coefficients N^ d and unf , while the parameters a q , (3 q , j q , 
6 q and Mh were taken to be flavour independent. For simplicity we also assumed that (kjjx = (kj_). With these 
choices, we were left with a total of 9 free parameters for the SIDIS-1 and SIDIS-2 fit parameterizations: 

iVj, Nj, Ng,, Ng { , a, p, 7 , S, M h . (18) 

Both fits gave good results. However, a study of the statistical uncertainties of the best fit parameters, and a 
comparison of the two sets of parameterizations, SIDIS-1 and SIDIS-2, clearly shows that SIDIS data are not presently 
able to constrain the large x behaviour of the quark (u, d) transversity distributions, leaving a large uncertainty in 
the possible values of the parameter /3. In fact, the range of Bjorken x values currently explored by HERMES and 
COMPASS experiments is limited to xb ^5 0.3. 

This uncertainty in the knowledge of the transversity distribution at large x values has relevant consequences 
when one uses the parameterizations extracted from SIDIS and e + e~ data for the study of single spin asymmetries 
in hadronic collisions. In this case the largest pion asymmetries have been measured at large Feynman x values, 
xp > 0.3; then, kinematical cuts imply that the transversity distribution is probed at even larger x values. 

In order to assess the possible relevance of the Collins effect in explaining a large value of An in pp collisions we 
should explore in greater details the large x contribution of the transversity distribution. We have then devised a 
simple analysis, to which we will refer to as the "scan procedure" and which is based on the following considerations. 
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• In our parameterizations, the large x behaviour of the quark tranvcrsity distributions is driven by the parameters 
f3 q , the exponents of the (1 — x) factor in M q {x) : see Eqs. (fTTj) and (|13p . Not surprinsingly, for the dominant u 
and d contributions, the values of the single f3 u — [3d — [3 parameter arc indeed very different in the SIDIS-1 and 
SIDIS-2 sets, despite the fact that they offer comparably good fits of SIDIS and e + e~ data. As a consequence, 
the two sets give strongly different estimates of the pion SSAs at large if in hadronic collisions. It is then 
natural to conclude that the choice of a flavour-independent /3 parameter, good for SIDIS data, is a much too 
strong assumption for the hadronic collisions, and must be released in this analysis. We also notice that both 
in SIDIS-1 and SIDIS-2 fits the best-fit values of <5f av = <5 un f = 5 are very close to zero. 

• We therefore start the scan procedure by performing a preliminary 9-parameter "reference fit" to SIDIS and e + e~ 
data taking, w.r.t. Eq. (fl8]l. [3 U ^ (3d and 8 — 0. We then let the two parameters (3 U and {3d vary independently 
over the range 0.0 — 4.0 by discrete steps of 0.5. Larger values of (3 would give negligible contribution to An- 
For each of the 81 points in this two-dimensional grid in /3d) space we perform a new 7-parameter fit to 
SIDIS and e + e~ data, keeping (3 Uj d fixed and 6 = 0, but leaving all remaining 7 parameters in Eq. (|18p free. 

• As a next step we select only those sets of parameters from the scan procedure over the {(3 Ul (3d) grid leading 
to an increment of the total x 2 of the fit, as compared to the corresponding 9-parameter reference fit, smaller 
than a given chosen value, A^ 2 . Notice that, since the reference fit and the scan fits have a different number of 
free parameters, the selection criterium is applied to the total \ 2 rather than to the \ 2 P er degree of freedom, 
Xdof ■ The chosen value of A% 2 is the same as that used in Refs. [HI, [lj| to generate the error band, following 
the procedure described in Appendix A of Ref. Q. It is worth noticing that all 81 points of our grid in {(3 U , (3d) 
lead to acceptable fits; this further confirms the observation that the SIDIS data cannot constrain the large x 
behaviour of the transversity distribution. 

• For each of the selected sets, we calculate the Collins pion SSA for polarized pp collisions, Eqs. (d])-©, in the 
kinematical regions of the available data from the E704 Collaboration at Fermilab and the STAR (for ir°) and 
BRAHMS (for charged pions) Collaborations at RHIC. 

Finally we generate a "scan band", by taking the envelope of ALL curves for An (it) obtained by using the sets 
selected in the scan procedure, and compare this band with the experimental data available. This band shows 
the potentiality of the Collins effect alone to account for An(p^p — > ttX) data while preserving a combined fair 
description (quantified by A^ 2 ) of the SIDIS and e + e~ data on Collins azimuthal asymmetries. 

• Notice that in the data sets used for the fits of the scan procedure we have also included the recent preliminary 
data by the COMPASS Collaboration on SIDIS off a transversely polarized proton target [76| . which were not 
available at the time when Refs. [l5l[l6j were published. Strictly speaking, therefore, the parameterizations used 
as starting point of the scan procedure are not the same as those published in Refs. fisl Il6|. However, we have 
verified that the new parameterizations are only slightly different from the previous ones and are qualitatively 
consistent with them. Therefore, here we will not present and discuss them further, referring for more details to 
a future complete upgrade of our parameterizations. Apart from the insertion of this new set of COMPASS data, 
all technical aspects of the fitting procedure followed here are the same as in Refs. mill, with the difference 
that, for the QCD evolution of the Collins function, limited to its collinear z-dependent D h / q (z) factor, we have 
attempted, in addition to an unpolarized-like evolution, also a transversity-like one, driven by the transversity 
evolution kernel. 



Results and comments 

We have computed An(p^p — > ttX) adopting, as explained above, a single set of collinear parton distributions [zl, 
[zH , two different sets for the pion collinear FFs [lH, [73| and two different (partial) evolution schemes for the Collins 
function. In Figs. 1-3 we show some of our results, avoiding the explicit presentation of other cases with very similar 
outcomes (some further comments are given below). For all results presented the Kretzer set for the unpolarizcd FFs 
and the unpolarized-like Collins evolution have been used. 

In Fig. Q] the scan band for An, as a function of xp at fixed scattering angles, is shown for charged pions and 
BRAHMS kinematics, while in Fig. [2] the same result is given, at fixed rapidity values, for neutral pions and STAR 
kinematics; analogous results, as a function of Pt at several fixed xf values, are shown for STAR kinematics in Fig. [3] 

These results allow to draw some first qualitative conclusions: 

• The Collins contribution to An is not as tiny as claimed in Ref. [l[; 

• The Collins effect alone might in principle be able to explain the BRAHMS charged pion results on An in the 
full kinematical range so far explored; 
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FIG. 1: Scan band (i.e. the envelope of possible values) for the Collins contribution to the charged pion single spin asymmetries 
An, as a function of xf at two different scattering angles, compared with the corresponding BRAHMS experimental data [35( |. 
The shaded band is generated, adopting the GRV98 and GRSV2000 sets of collinear PDFs, the Kretzer FF set and an 
"unpolarized-like" evolution for the Collins function, following the procedure explained in the text. 
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FIG. 2: Scan band (i.e. the envelope of possible values) for the Collins contribution to the neutral pion single spin as ym metry 
An, as a function of xf at two different rapidity values, compared with the corresponding STAR experimental data [3J]. The 
shaded band is generated, adopting the GRV98 and GRSV2000 sets of collinear PDFs, the Kretzer FF set and an "unpolarized- 
like" evolution for the Collins function, following the procedure explained in the text. 
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• The full amount of the ir° STAR data on An cannot be explained by the Collins contribution alone. The Collins 
effect might be sufficient for the small xf portion of the data; however, it is not sufficient for the medium-large 
xp range of STAR data, xf > 0.3. 

The results obtained with a different choice of the fragmentation functions (the DSS set) are qualitatively very 
similar in the large xf regions. They are instead smaller in size at smaller xf, due to the large gluon contribution in 
the leading order (LO) DSS fragmentation functions. The use of a transversity-like Collins evolution, rather than the 
unpolarizcd one, does not lead to any significant difference, in all cases. 

At this point, in order to fully assess the role of the Collins effect in understanding the large SSAs for neutral pions 
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FIG. 3: The same as for Fig. [2] but with the STAR data plotted vs. the pion transverse momentum, Pt, for different bins in 
x F , x F = 0.28, 0.37, 0.43 and 0.50. 



measured at large xf by the STAR Collaboration at RHIC, we have performed several further tests. 

First of all, we should make it clear that the scan bands presented in our plots have nothing to do with the statistical 
error bands presented in Refs. flH [l6j . There, the error bands are generated by estimating the uncertainty in the 
best fit values of the parameters, according to the procedure described in detail in Appendix A of Ref. Q. Instead, 
the scan bands in this paper are obtained by simply taking the envelope of all curves generated by the selected best 
fit sets within the full grid in f3 u ,d- 

It is not clear how to combine the statistical error band, associated with the full 7 or 9 free-parameter best fits of 
SIDIS and e + e~ data, with the scan bands. Therefore, in order to understand to what extent the statistical errors 
on the best fit parameters may affect the capability of the Collins effect to reproduce the large x f STAR data, we 
have adopted the following strategy: besides considering the envelope of the full set of curves produced by the scan 
procedure, we have considered explicitly each of these curves, isolating the set leading to the largest asymmetries in 
the large xf region; we have then evaluated, as in Appendix A of Ref. Q, the corresponding statistical error band, 
which covers larger values of the asymmetry. Our result is presented in Fig. [4] Again, it appears that the Collins 
effect alone cannot account for the large xf data. Notice also that trying to fit the large xf data on An might lead 
to an overestimation of the same data at smaller xf, which have tinier error bars. 
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FIG. 4: The Collins contribution to the neutral pion single spin asymmetry Am, compared with the corresponding STAR 
experimental data at two fixed pion rapidities [34|]. The shaded statistical error band is generated, adopting the GRV98 and 
GRSV2000 sets of collinear PDFs, the Kretzer FF set and an "unpolarized-like" evolution for the Collins function, starting 
from the 7-parameter fit in the grid procedure that maximizes the neutral pion SSA in the large if region and applying the 
error estimate procedure described in Appendix A of Ref. ||. See text for more details. 



There is still another issue that deserves some attention. Although simplified, our parameterization of the TMDs 
and of their functional shape involves in principle a huge number of fit parameters. Since most of these parameters 
are highly correlated, adopting larger set of parameters would lead to larger uncertainties in their value. Therefore, 
reasonable fits require a reduction in the number of fit parameters and involve a careful choice of the most significant 
ones. This choice may have consequences on the allowed values of the asymmetries, particularly for kinematical regions 
not covered by the data sets used for the fitting procedure. In the present scan procedure we use a 7-paramctcr fit 
and a grid of values for the two additional parameters /3 U , fid- In order to investigate if a larger set of free parameters 
for the scan procedure could modify our conclusions about the Collins effect for the STAR data in the large xf range, 
we have repeated our scan procedure by starting from a preliminary reference fit with 13 free parameters, 

a u , a d , f3 u , (3 d , N fav , N un{ , 

Tfav: Tunf; ^fav; ^unf; 

M h , (19) 

and generating again the scan band on the two dimensional grid for the fixed f3 u _d parameters by fitting for each 
grid point the remaining 11 parameters. This naturally results in a sizably larger scan band, with its upper edge 
approaching better An at the larger xf values. However, also in this case, looking at all the 81 fit sets we find that 
the curve with the best behaviour at large xf approaches the upper edge of the scan band in the full xf range. It 
therefore largely misses (overestimates) the lower x f values of the asymmetry. In Fig. [5] we present the statistical 
error band on the Collins contribution to A^{p^p — > n°X) generated, following Appendix A of Ref. Q, from the 
11-paramctcr best- fit set which optimizes the agreement with STAR data at large xp- 

Let us finally make some comments on the (charged and neutral) pion SSAs for the E704 kinematics [HI, H^] . The 
situation in this case is complicated by the fact that, contrary to the STAR and BRAHMS kinematics, the unpolarized 
cross sections are largely underestimated within the TMD LO factorized approach when adopting the values of (k^_) 
and (p 2 ^ extracted from SIDIS data. Indeed, much larger effective values are required to reconcile the TMD estimates 
with data, as it was shown in Ref. [ijj. However, this fact should have less influence on the SSAs, defined as ratios 
of (sums and differences of) single-polarized cross sections. 

We have therefore directly applied the scan procedure, as illustrated above, also to the E704 results. Again, it turns 
out that the scan band could cover the data for the neutral pion SSA, with some problem at the largest xf values. 
However, it largely misses the huge charged pion SSAs observed in the same kinematical region. 
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FIG. 5: The same as for Fig. [4] but this time with the statistical error band generated starting from the scan procedure with 
11 free parameters. See text for more details. 

III. CONCLUSIONS 

We have investigated the possible role of the Collins effect in explaining the large SSAs observed in p^p — > ttX 
reactions; we have done so within a TMD factorized scheme, and have revisited a previous work on the same is- 
sue, correcting a numerical error and using new experimental data and new phenomenological information on the 
transversity distribution and the Collins function. 

We can conclude that, to the best of our present knowledge, based on SIDIS and e + e~ data, the Collins effect alone 
seems to be able to reproduce the available RHIC data on pion single spin asymmetries in polarized pp collisions, 
only in the small Feynman x region, xf < 0.3. Above that, which is the region where the values of An increase, the 
Collins effect alone is not sufficient. 

Additional mechanisms are required in order to explain the size of the Ajq asymmetry in this region. One can 
obviously think of the Sivers effect [H, US SH HH ■ Since TMD factorization has not been proven and is still under 
active debate for single inclusive particle production in hadronic collisions and since universality breaking effects are 
possible, one does not know exactly how to use the parameterizations of the quark Sivers functions extracted from 
SfDIS data in pp collisions. A recent use of the SIDIS Sivers functions in a collincar higher-twist approach to SSAs 
in pp collisions - rather than in the TMD factorized approach - has been found to give a sizable contribution to Ajv, 
but with the wrong sign (77J ■ A similar conclusion holds in a modified generalized parton model approach with TMD 
factorization [78[ . Notice that such a problem does not occur if one simply adopts the SIDIS Sivers functions in the 
TMD factorized scheme [7§, H(| ■ Much further investigation is necessary. 
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